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ABSTRACT 

Despite the many studies on stellar nucleosynthesis published so far, the scenario for the 
production of Cu in stars remains elusive. In particular, it is still debated whether copper 
originates mostly in massive stars or type la supernovae. To answer this question, we compute 
self-consistent chemical evolution models taking into account the results of updated stellar 
nucleosynthesis. By contrasting copper evolution in ui Cen and the Milky Way, we end up 
with a picture where massive stars are the major responsible for the production of Cu in 
uj Cen as well as the Galactic disc. 

Key words: Galaxy: evolution - globular clusters: individual: ui Centauri - nuclear reactions, 
nucleosynthesis, abundances - stars: abundances. 



1 INTRODUCTION 

Almost two decades ago, Sneden & Crocker (1988) and Sneden, 
Gratton & Crocker (1991) discussed the possible mechanisms for 
the synthesis of Cu in stars in the light of their high-resolution ob- 
servations (showing a decline in [Cu/Fe] with decreasing [Fe/H] in 
Galactic stars) and favoured a secondary origin for Cu through the 
weak component of the .s-process in massive stars. 

The .s-process has long been thought to contribute to the pro- 
duction of Cu in stars (see, e.g., Burbidge et al. 1957). The main 
i-process, operating in low-mass asymptotic giant branch (AGB) 
stars, does not contribute more than 5 per cent to the synthesis of Cu 
isotopes (Raiteri et al. 1993; Travaglio et al. 2004) and will be ne- 
glected hereinafter. The weak .s-process in massive stars depends on 
metallicity. This mechanism, in fact, requires neutron fluxes origi- 
nating mainly from the reaction 22 Ne (a, n) 25 Mg, where the large 
abundance of 22 Ne during core He burning derives from the orig- 
inal CNO nuclei trasmuted into 14 N in the H burning ashes, fol- 
lowed by double a-capture on 14 N. The residual 22 Ne, left behind 
at core He exhaustion, is fully consumed in the O-rich zone during 
convective shell C burning. A small primary yield of Cu, 5 to 10 per 
cent of its solar abundance depending on the location of the mass 
cut, derives from explosive nucleosynthesis in the inner regions of 
core-collapse supernovae (SNe II). 

In the early nineties, only 25 per cent of solar Cu was at- 
tributed to the weak ^-process in massive stars (Raiteri et al. 1993). 
It was advanced then (Matteucci et al. 1993; Mishenina et al. 2002) 
that most solar Cu originated in type la supernovae (SNe la). Cur- 
rent SNIa models, however, predict a negligible production of Cu 
during thermonuclear explosions (Iwamoto et al. 1999; Travaglio, 
Hillebrandt & Reinecke 2005). 

It has been recently claimed that up-to-date, high-resolution 
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spectroscopic data for solar neighbourhood stars belonging to dif- 
ferent Galactic substructures (halo, thick-disc and thin-disc), bulge- 
like stars, globular clusters (GCs) and dwarf spheroidal galaxies 
(dSphs) can all be well understood if Cu is mostly a secondary 
product of massive stars, with a small primary contribution from 
explosive nucleosynthesis (Bisterzo et al. 2005). While this con- 
clusion simply rests on the comparison of extant observations with 
present stellar nucleosynthesis (and Cu yields from stars ultimately 
evolving to SNe II must still be viewed with great care - see, e.g., 
Heger et al. 2001; Hoffman, Woosley & Weaver 2001; Rauscher 
et al. 2002), a detailed description of copper evolution in different 
galactic structures is still missing from the literature. 

In this Letter, we try to fill in this gap by analysing the evolu- 
tion of Cu in the framework of detailed galactic chemical evolution 
(GCE) models in which Cu originates mostly from (i) massive stars 
or (ii) SNe la. In particular, we contrast the evolution of copper in 
the Milky Way with that in the anomalous GC u Cen. We conclude 
that scenarios where copper is produced primarily in SNe la are 
no longer supported by the observations. This imposes strong con- 
straints on the astrophysical processes responsible for the synthesis 
of Cu in stars. 



2 MODELS VERSUS OBSERVATIONS: IMPLICATIONS 
FOR STELLAR YIELDS AND CHEMICAL 
EVOLUTION TIME SCALES 

Nowadays, high-resolution spectroscopic data of [Cu/Fe] are avail- 
able for large samples of field Galactic stars (see Fig. 1, left panel). 
GCs generally follow the trends defined by halo field giants (Sim- 
merer et al. 2003): a flat distribution, [Cu/Fe] = -0.75 ± 0.2, 
for low-metallicity stars up to [Fe/H] ~ —1.8 dex, followed by 
a linear increase with a slope close to 1 in the metallicity range 
-1.5 < [Fe/H] < -1 (GC data are not plotted in Fig. 1). At 
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Figure 1. Left panel: [Cu/Fe] in Galactic disc and halo stars, as measured by Prochaska et al. (2000; big empty squares), Mishenina et al. (2002; stars), Reddy 
et al. (2003; small filled squares), Allende Prieto et al. (2004; small rilled upside-down triangles), Bihain et al. (2004; big empty triangles) and Reddy et al. 
(2006; dots) in the solar neighbourhood. Also shown are [Cu/Fe] values for the very metal-poor stars BD+17° 3248 (Cowan et al. 2002; big filled square), 
HD 1 15444, HD 122563 (Westin et al. 2000; big filled triangles) and the extremely metal-poor giant CS 22892-052 (Sneden et al. 2003; big star) and TP-AGB 
star CS 30322—023 (Masseron et al. 2006; pentagon). Right panel: [Cu/Fe] in the globular cluster w Cen. Data from Cunha et al. (2002) and Pancino et al. 
(2002; filled circles) have been homogenized as discussed in Romano et al. (2007). Superimposed on the data are our model predictions for uj Cen (right 
panel, thin lines) and for the solar vicinity (left panel, thick lines), in cases where Cu originates mostly from SN explosions (Models 1 and 2, Table 1 , long- 
and short-dashed lines, respectively; both type fl and type la SN explosions are considered). The effect of completely removing the (primary) explosive SNfl 
contribution is also shown as a dotted line (only in a model for the solar vicinity; left panel, Model 4, Table 1). See text for details on different models' 
prescriptions. 



[Fe/H] > —0.8, [Cu/Fe] jumps above the solar value. Then, a 
'bending' appears for disc stars with —0.8 < [Fe/H] < (Reddy 
et al. 2003), though there is a hint that [Cu/Fe] might start increas- 
ing again at higher metallicities (Allende Prieto et al. 2004). 

The Galactic GC u) Cen stands as a notable exception. The 
[Cu/Fe] ratios of its most metal-rich stars are definitely lower than 
the Galactic trend (Cunha et al. 2002; Pancino et al. 2002; see also 
Figs. 1 and 2), which can be understood as a shift in the [Cu/Fe] 
relation to higher [Fe/H], The only other systems known to have 
unusually low [Cu/Fe] ratios are the Sagittarius dSph (McWilliam 
& Smecker-Hane 2005) and the Large Magellanic Cloud (LMC; 
Pompeia, Hill & Spite 2005). In particular, the similarity between 
the Cu values of uj Cen and Sagittarius adds to the list of common 
chemical peculiarities (spread in [Fe/H] values, strong enhance- 
ment of i-process elements) which suggest that the two systems 
followed similar chemical enrichment histories, thus supporting the 
idea that uj Cen is the surviving nucleus of an accreted dwarf galaxy 
(McWilliam & Smecker-Hane 2005). 

Superposed on the data in Figs. 1 and 2 are different model 
predictions for uj Cen (thin lines) and the Milky Way (thick lines). 
These latter refer to the solar neighbourhood region which, in the 
framework of the adopted model, is represented by a ring 2 kpc 
wide centred in the Sun. The model of chemical evolution for the 
Galaxy is the two-infall model of Chiappini, Matteucci & Gratton 
(1997) and Chiappini, Matteucci & Romano (2001) - where details 



about model assumptions and basic equations can be found, except 
for the adopted stellar lifetimes (which are now taken from Schaller 
et al. 1992) and stellar yields (see the detailed discussion below). 
The model of chemical evolution for uj Cen is the one described in 
Romano et al. (2007), assuming that wCen is the remnant of an 
ancient nucleated dSph evolved in isolation and then accreted by 
the Milky Way. We adopt this simple model even if a recent study 
by Villanova et al. (2007) suggests that the chemical evolution of 
uj Cen might be more complex, with a not unique age-metallicity 
relation. In fact, accounting for the observed spread needs a full 
dynamical treatment, while the main conclusions on Cu nucleosyn- 
thesis in stars are likely to remain the same. 

As regards the adopted stellar nucleosynthesis: 

(i) the s-process yields in massive stars vary as functions of both 
the initial mass and metallicity of the stars; we follow the prescrip- 
tions of either Matteucci et al. (1993; Models 1, 2 and 4, Table 1; 
long-dashed, short-dashed and dotted lines in Fig. 1, respectively) 
or Nomoto et al. (2007; Model 3, Table 1; solid lines in Fig. 2); 

(ii) explosive yields of Cu of individual core-collapse SNe 
are taken from either Matteucci et al. (1993, their model M, 
M^uSNcii = 7 x 10~ 6 Mq, independent of mass and metallicity; 
Models 1 and 2, Table 1; long- and short-dashed lines in Fig. 1, 
respectively) or Nomoto et al. (2007; Model 3, Table 1; solid lines 
in Fig. 2); 

(iii) the amount of Cu in the ejecta of SNe la is of 2 x 10 -4 M© 
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Figure 2. Same data as Fig. 1. The [Cu/Fe] versus [Fe/H] relation in us Cen is now compared to that for solar neighbourhood stars. Model predictions for us Cen 
(thin solid line) and the solar vicinity (thick solid line) now refer to the evolutive picture where most of solar Cu comes from the weak .v-process in massive 
stars (stellar nucleosynthesis prescriptions as in Model 3, Table 1). See text for details. 



Table 1. Prescriptions on stellar Cu production adopted by different GCE models. 



Model 


M cj 

Cu weak 




M cj 

m Cu SNoII 


M cj 

Cu SNoIa 


Notes 


1 


Matteucci et al 


1993 


Matteucci et al. 1993 (their model M) 


Matteucci et al. 1993 (their model M) 


Long-dashed lines in Fig. 1 


2 


Matteucci et al 


1993 


Matteucci et al. 1993 (their model M) 


Iwamoto et al. 1999 (their model W7) 


Short-dashed lines in Fig. 1 


3 


Nomoto et al. 


2007 


Nomoto et al. 2007 


Iwamoto et al. 1999 (their model W7) 


Continuous lines in Fig. 2 


4 


Matteucci et al 


1993 





Matteucci et al. 1993 (their model M) 


Dotted line in Fig. 1 



(Matteucci et al. 1993, their model M; Models 1 and 4, Table 1; 
long-dashed and dotted lines in Fig. 1) or two orders of magnitude 
lower (Iwamoto et al. 1999, their model W7; Models 2 and 3, Ta- 
ble 1 ; short-dashed lines in Fig. 1 and solid lines in Fig. 2), and it 
is independent of metallicity. 

Total yields of copper from massive stars by Matteucci et al. (1993, 
squares) and Nomoto et al. (2007, dots), comprising both the weak 
component and the explosive one, are shown and compared one 
another in Fig. 3, for solar metallicity stars (upper panel) and zero 
metallicity objects (lower panel). At Z = 0, the yields reflect the 
pure explosive contribution. It is seen that, while the explosive 
yields adopted by Matteucci et al. (1993) do not vary with the ini- 
tial mass of the star, those computed by Nomoto et al. (2007) vary 
steeply as a function of mass, especially for Mini > 20 Mq. 

Models where most of solar Cu comes from SNIa explo- 
sions following the prescriptions of Matteucci et al. (1993) match 
very well the observational data for solar neighbourhood stars (see 
Fig. 1, left panel, thick long-dashed line), but leave the lower 
[Cu/Fe] ratios in uj Cen totally unexplained (see Fig. 1, right panel, 
thin long-dashed line). When reducing the contribution from SNela 
according to the results of Iwamoto et al. (1999), a better agreement 
with the relation observed for us Cen members is found (except, per- 



haps, for the two objects lying at the highest metallicities - Fig. 1, 
right panel, thin short-dashed line), but the agreement with the so- 
lar neighbourhood data at [Fe/H] > —1.0 is completely destroyed 
(Fig. 1, left panel, thick short-dashed line). If, instead, most of solar 
Cu is produced through the weak s-process in massive stars, both 
the solar neighbourhood and us Cen relations are satisfactorily re- 
produced (Fig. 2, thick and thin solid lines). It is worth stressing 
at this point that, in order to best fit the solar abundance of Cu in 
the framework of our GCE model, we had to slightly increase (by 
a factor of 1.5) the yields of Cu by Nomoto et al. (2007) for so- 
lar metallicity stars. The applied correction factor, however, is well 
inside current uncertainties on the weak .s-process modelling (see 
Rauscher et al. 2002 for a discussion on this subject). It is also 
worth emphasizing how dramatically important the production of 
Cu from core-collapse SNe is in the very early evolution of the 
Galactic halo. The SNII contribution sets the level of the 'plateau' 
observed at [Fe/H] < —1.8. This can be clearly seen from Fig. 1 
(left panel), where the thick dotted line (Model 4, Table 1) shows 
the effect of removing any primary Cu production from explosive 
nucleosynthesis in massive stars from the Milky Way model. In this 
case, the 'plateau' at low [Fe/H] is no longer expected, at variance 
with the observations, and only the steep rise towards high [Cu/Fe] 
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Figure 3. Stellar yields of Cu (in units of solar mass) for solar-metallicity 
and zero-metallicity massive stars from Matteucci et al. (1993, squares) and 
Nomoto et al. (2007, dots) as functions of the initial mass of the stars. At 
solar metallicity (upper panel), the yields comprise a weak s-process com- 
ponent as well as an explosive one, while for Z = (lower panel) only the 
explosive component is present. 



at high metallicities can be reproduced. Furthermore, notice that 
when adopting the yields of Nomoto et al. (2007) rather than those 
of Matteucci et al. (1993), a lower plateau level is obtained (Fig. 2), 
because of the lower amount of Cu in the ejecta of low-metallicity 
SNell predicted by those authors (see Fig. 3, lower panel). 
Now we can draw the following conclusions: 

(i) After a short phase in which the primary contribution from 
explosive nucleosynthesis in core-collapse SN dominates, the evo- 
lution of Cu in galaxies of different type is regulated mostly by the 
weak i-process occurring in massive stars. 

(ii) Up-to-date, metallicity-dependent stellar yields of Cu by 
Nomoto et al. (2007) including both a primary contribution from 
explosive nucleosynthesis and a secondary contribution from the 
i-process in massive stars allow us to reproduce very well the ob- 
served behaviour of [Cu/Fe] versus [Fe/H] in the solar vicinity as 
well as in uj Cen. 

(iii) When restricting our analysis to a single environment (the 
solar vicinity) a degeneracy is found in the solution, in the sense 
that the decline in [Cu/Fe] with decreasing [Fe/H] can be ascribed 
to either the reduced extent of the weak s-process in massive stars at 
low [Fe/H], or the delayed Cu production from SNela. Only study- 
ing objects which have chemically enriched along different evolu- 
tionary paths allows us to break the degeneracy. 

Finally, in Fig. 4 we compare the abundances of Cu measured 
in a sample of bulge-like stars by Pompeia (2003) with our model 
predictions for the inner Galaxy (bulge-like stars are thought to 
originate near the Galactic bulge). A substantial contribution to the 
Galactic Cu production from SNela (nucleosynthesis prescriptions 



Figure 4. Predictions from Models 1 (long-dashed line), 2 (short-dashed 
line) and 3 (solid line) computed for the inner Galaxy (Rq = 4 kpc) com- 
pared to measurements of [Cu/Fe] in a sample of bulge-like stars (data from 
Pompeia 2003). 

as in Model 1, Table 1 ; long-dashed line) is clearly ruled out, while 
models with a minimum Cu production from SNela can well fit the 
data (Models 2 and 3, short-dashed and solid lines, respectively). 
We conclude that the study of Cu abundances in different environ- 
ments, from both an observational and a theoretical point of view, 
is of fundamental importance in order to understand the origin of 
Cu in the universe. 
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